Time Evolution of Entropy of a Charged Domain Wall during Gravitational Collapse 
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We study the time evolution of the entropy of a collapsing charged spherical domain wall, from 
the point of view of an asymptotic observer, by investigating the entropy of the entire system 
(i.e. charged domain wall and radiation) and induced radiation alone during the collapse. By taking 
the difference, we find the entropy of the collapsing charged domain wall, since this is the object 
which will form a black hole. We find that for large values of time (times larger than tf-/R+ » 8), 
the entropy of the collapsing domain wall is a constant, which is of the same order as the Bekenstein- 
Hawking entropy. 
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I. INTRODUCTION 

Hawking's 1974 paper (see Refs. [H, [1|) sparked great 
interest in both the existence of the radiation given off 
during the loss of mass for a pre-existing static black 
hole and speculations as to what his result implied for 
the validity of quantum field theory in the presence of a 
black hole. Hawking showed that the radiation given off 
was in fact thermal, which only increased the interest. 

In 1972 Bekenstein argued that a black hole of mass 
M has an entropy proportional to its area (see Ref. 0). 
Further calculations by Gibbons and Hawking showed 
that the entropy of a black hole is always a constant, re- 
gardless of the type of metric which is used (see Ref. d|). 
They showed that the expression for the entropy is given 
by 



~>BH 



(1) 



where A] lor is the surface area of the horizon. 

Recently, subsequent work has been done in theories 
of quantum gravity (such as String Theory and Quan- 
tum Loop Gravity) to reproduce this result (see for ex- 
ample Refs. 0, So far, most of these calculations 
have involved only pre-existing black holes. The typical 
method for finding the entropy of the black hole is to 
first determine the temperature of the black hole using 
the Bogolyubov method. Here, one considers that the 
system starts in an asymptotically flat metric, then the 
system evolves to a new asymptotically flat metric. By 
matching the coefficients between the two asymptotically 
flat spaces at the beginning and end of the gravitational 
collapse, the mismatch of these two vacua gives the num- 
ber of particles produced during the collapse. Hence, 
what happens in between is beyond the scope of the Bo- 
golyubov method. Therefore the time-evolution of the 
thermodynamic properties of the collapse cannot be in- 
vestigated in the context of the Bogolyubov method. 

In this paper we will investigate the collapse of a 
charged spherically symmetric infinitely thin domain wall 
(representing a shell of charged matter). The purpose 
of this paper is investigate the time-dependent entropy 
evolution of the charged domain wall during gravitational 
collapse and verify the claim that the entropy of the black 
hole is independent of the background metric. There- 



fore we will compare the results to Eq. ([T]). The details 
about the collapse will depend on the particular folia- 
tion of space-time used to study the system. For the 
purposes of the question we are studying, we are out- 
side asymptotic observers. Thus, we study the collapse 
of the domain wall from the view point of a stationary 
asymptotic observer. 

In order to study the time-dependent thermodynamic 
processes of the system, we will repeat the steps first de- 
veloped in Ref. [7| for the massive case. This is done by 
first developing the partition function, in Section IIV Al 
To do so we shall employ the so-called Liouville-von Neu- 
mann approach (see Refs. [1, 0). This approach utilizes 
the invariant operator approach developed by Lewis and 
Reisenfeld (see Ref. [l(|). In the invariant operator for- 
malism, it was observed that any Hcrmitian operator I(t) 
satisfying the Liouville-von Neumann equation, 
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can be used to find the exact wavefunctions of the 
Schrodinger equation, up to an overall time-dependent 
phase factor. In the Liouville-von Neumann approach, 
the time-dependent Hamiltonian is replaced by the time- 
dependent invariant operator, since the eigenvalues of 
this operator are time-independent and all the time- 
dependence of the system is placed into a non-linear aux- 
iliary equation. 

To find the parameter j3, we shall use the so- 
called Functional Schrodinger formalism (sec for exam- 
ple Refs. [lUil), in Section HHCl In general, the Func- 
tional Schrodinger formalism yields a functional differ- 
ential equation for the wavefunctional, ^[g^, X 11 , O], 
where g^ is the metric, X^ the position of the object, 
$ is a scalar field and O denotes all the observer's de- 
grees of freedom. Since the Functional Schrodinger for- 
malism depends on the observer's degrees of freedom, 
one can introduce the "observer" time into the quan- 
tum mechanical processes in the form of the Schrodinger 
equation. The wavefunctional is then dependent on the 
chosen observer's time, hence one can view the quantum 
mechanical processes of a given system under any choice 
of space-time foliation. One benefit of this formalism 
is that one can solve the time-dependent wavefunctional 
equation exactly using the invariant operator method, in 
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particular in the present paper. Here we study the prop- 
agation of induced radiation, represented as a scalar field 
in the background of gravitational collapse, which is gov- 
erned by the harmonic oscillator equation with a time- 
dependent frequency. We solve the equations of motion 
to find the time-dependent wavefunctional for the sys- 
tem. One can then define the occupation number of the 
induced radiation using the wavefunctional, which is the 
gaussian overlap between the initial vacuum state and the 
wavefunction at any given later time. By taking the oc- 
cupation number to be of the form of the Planck distribu- 
tion, we can then define (3 which is also time-dependent. 

Since the partition function is time-dependent, one 
can then study the time-dependent thermodynamic pro- 
cesses associated with gravitational collapse. Therefore 
the Functional Schrodingcr formalism and Liouville-von 
Neumann approach together go beyond the approxima- 
tions of the Bogolyubov method, since the system is al- 
lowed to evolve over time, which allows one to investigate 
the intermediate regime during the collapse. 



II. MODEL FOR COLLAPSING CHARGED 
DOMAIN WALL 

To study a realization of the formation of a black hole, 
we consider a spherical Nambu-Goto domain wall (repre- 
senting a domain wall of matter) that is collapsing. This 
will be done using the Functional Schrbdinger equation, 
where we will consider a minisupcrspacc version of the 
Wheeler-de Witt equation. Since we are interested in 
the entropy of the collapsing domain wall from the point 
of view of an asymptotic observer, the relevant degree of 
freedom for the collapsing spherical domain wall is the 
radial degree of freedom R(t). The metric of the system 
is then chosen to be the solution to Einstein's equations 
for a spherical domain wall. Therefore, outside the wall 
the metric is Schwarzschild, as follows from the spherical 
symmetry 



ds 2 
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where 



dtt 2 = dB 2 + sin' 



(3) 



(4) 



By Birkhoff 's theorem, the metric inside the domain wall 
must be flat, hence Minkowski 



ds 2 



-dT 2 + dr 2 + r 2 dn 2 



(5) 



where T is the interior time. The interior time is related 
to the asymptotic observer time t via the proper time r 
of the domain wall. 
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and 



where 




2GM Q 2 



R 



R 2 ' 



By taking the ratio of these equations we then have 



(7) 



(8) 



(9) 



where R = dR/dt. 

By integrating the equation of motion for a spherically 
symmetric domain wall, Wang et. al. and Lopez (see 
Refs. [THHH) found that the mass is a constant of motion 
and is given by 

Q 2 




M = AnaR 2 



t/1 + RJ - 2naGR 



2R 



(10) 



where a is the surface tension of the domain wall and 
R T = dR/dr. Before we proceed we wish to discuss the 
physical relevance of Eq. (fTU|). First consider the case 
where R T = 0, i.e. for a static domain wall. The first 
term in the square bracket is just the total rest mass 
of the domain wall. When the domain wall is moving, 
i.e. R T ^ 0, the first term in the square bracket takes 
the kinetic energy of the domain wall into account. The 
last term in the square bracket is the self-gravity, or the 
binding energy of the domain wall. Finally, the last term 
in Eq. (TTU)) is the electromagnetic contribution to the 
total mass (energy). Therefore we can identify Eq. (p~0|) 
as the total energy of the system, hence the Hamiltonian 
of the system. Thus, we will refer to Eq. (|10l) as the 
Hamiltonian. 

Using Eq. (|Tj) the Hamiltonian for the wall, Eq. (fTOl) . 
can be written as (see Ref. [TJ[ ) 



H = 47TCr/ 3/2 i? 2 



1 



- 2-kuGR 



p-(l-f)R 2 Jf 2 -R 2 



(11) 

where R = dR/dt. The canonical momentum near the 
horizon, the R ~ Rh regime, is given by 



4tthR 2 R 



V7\// 2 - R 2 



where 



fx=l- 2ixoGRh 



8tt<tR 3 h ' 



(12) 



(13) 



In this region the wall Hamiltonian Eq. (TTTT) in terms of 
the canonical momentum is then 



^lipl 2 R 2 



' p-R 2 

V(/n i? ) 2 +/(47T M i?2)2 (14) 
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The Hamiltonian has the form of the energy of a rel- 
ativistic particle, \/p 2 + m 2 , with a position-dependent 
mass. 

In the region R ~ R H , Eq. © tells us that / ~ 0. 
From Eq. (fl"4"]l one can see that the mass term can be 
neglected, thus Eq. (fT4|) reduces to 



(15) 



where we have chosen the negative sign since the domain 
wall is collapsing. 

Considering only the classical solution, the near- 
horizon solution of Eq. can be written as 



R. 



-f 



1 



2GM Q 2 



R R 2 

For the nonextremal case, we consider the equality 
2GM Q 2 



1 



R H R 2 



0. 



where Rh is given by 



R H = GM ± y/(GM) 2 - Q 2 . 



(16) 



(17) 



(18) 



The plus sign is the outer and the minus sign is the inner 
horizon. To distinguish them we write 



R+ = GM + y/(GM) 2 - 



R- = GM - y/{GM) 2 - Q 2 
Therefore we can write Eq. (jH]) as 



(19) 



(20) 



R 



/= 1-^ 1-^ 



R 



R. 



R 



/+/- 
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Since in the near horizon limit / — > 0, we can work with 
two different limits, either / + — ► or /_ — > 0. In both 
cases R w ±/ + /_ — > 0. For an asymptotic observer 
watching the collapse, the horizon of interest is / + . For 
/ + — > 0, /_ goes to a finite constant number, thus we 
have that the solution for Eq. (|16p is then 



R(t) w i?+ + (i? - i?+)e 



-f-t/R+ 



(22) 



Eq. implies that a collapsing domain wall crosses its 
own horizon only after an infinite amount of asymptotic 
observer time t. 

Here we note that Eq. (|22[) implies that in the case 
Q 2 > (GM) 2 , the argument in the exponential becomes 
imaginary. Thus the solution for the collapsing charged 
domain wall contains sinusoidal terms which oscillate 
about a radius R outside of the horizon. In this case 
the horizon never forms, thus never creating a black hole. 
Therefore, we shall ignore this solution and only investi- 
gate solutions where (GM) 2 > Q 2 . 



III. INDUCED RADIATION 

Here we consider the radiation given off during gravi- 
tational collapse. We will consider the radiation given off 
by the entire system and the particles produced during 
collapse alone. 

To investigate the radiation, we consider a complex 
scalar field $ in the background of the collapsing domain 
wall. The complex scalar field is decomposed into a com- 
plete set of basis functions denoted by {/&(?*)} 



^2(a k (t)Mr)+ib k (t)f k (r)). (23) 



The exact form of the functions fk(r) will not be im- 
portant for us. We will, however, be interested in the 
wavefunction for the mode coefficients {a k } and {bk}- 

The Hamiltonian for the complex scalar field modes is 
found by inserting Eq. ((23]) and Eq. © and Eq. © into 
the action 



5$ = 



-g-cT (D^T 



where 



iqAf, 



(24) 



(25) 



is the covariant derivate, q is the charge of the field and 
is the vector potential. However, from Ref. [l5[, we 
can write the vector potential for the collapsing charged 
shell as 



(AO 



Q 



0,0,0 



(26) 



The Hamiltonian for the complex scalar field modes, ar- 
rived at from Eq. (|24[) . takes the form of uncoupled sim- 
ple harmonic oscillators for each mode with independent 
mass and couplings due to the non-trivial metric. Using 
the principle axis transformation, the Hamiltonian can 
be diagonalized and written in terms of eigenmodes de- 
noted by c and d, where c is the eigenmodc for a and d is 
the eigenmode for b. It was shown in Ref. [l6| that in the 
regime of R ~ the Hamiltonian for a single c mode 
can be written as 



H.-f 1 ^ | j qQy 
2m 2m 



and for a single d mode 



Hd = f K_ i 29v 

2m 2m 



K 



-d 2 



(27) 



(28) 



where n c is the momentum conjugate to c, is the 
conjugate momentum to d, and where m, y and K are 
constants whose precise values are not important to us. 

The physical meaning of the c and d modes are inferred 
from the Noether current. From Ref. [l6| . we can write 
the Noether current as 



dd) + qQy (c 2 + d 2 ) 



(29) 
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where c = dc/dt and d — dd/dt. The mode coefficients 
satisfy the commutation relations given by 



[c,c] = [d,d] = 1, 



(30) 



Thus we treat the background classically and the radia- 
tion quantum mechanically. 

Dividing Eq. (|36j) through by /, we now rewrite in the 
standard form of a harmonic oscillator 



hence we see that the the Noether current satisfies the 
commutation relations 



[iV<j,,c] = c, and [iV$,<i] 



(31) 



Therefore we see that the c modes create particles with 
+q worth of charge, while the d modes create particles 
with — q worth of charge. 



A. Entire System 



From Eq. (|T5J) and Eqs. <J27J) and ([28]) we can write the 
Hamiltonian of the entire system as 



H — H wa u + Hb 



■ qQy 



,n 



.qQy , K 



f R + 2 m 2m ^ 2 ^ 2m 2m 2 



n 



K 



(32) 



where 11^ = —id/dR, TL r = —id/dc and II^ = —id/dd. 
As pointed out in Ref. [16|, the total wavefunction is in- 
dependent of the charge charge of the modes. The wave- 
function for the entire system is then a function of c, d, 
R, and t, which we can write as 



# = f(c, d,R,t). 



(33) 



Using the standard quantization procedure and 
Eq. (|32|). we can then write the Functional Schrodingcr 
equation as 



if 



/ a 2 * / d 2 ^ K 



dR 2m dc 2 2m dd 2 



— d 2 ^ = i—. 

2 dt 
(34) 

To solve Eq. (|34p we choose the classical background for 
the collapsing domain wall. Since the distance of the 
domain wall only depends on time, see Eq. (|16j) . we can 
then write 



dt f a 2 * /a 2 * a 2 

j f - t I c v[r -| d ^ = i 

1 dR 8t 2m dc 2 2m 3d 2 2 



K 



.a* 

~dt 



Subtracting the first time from both sides and grouping 
like terms, we can then write Eq. (|34[) as 



J_d^__J_cP^_ K 2 A 2 _ .0* / _ dt\ 
2m dc 2 2m dd 2 + 2 + 2 dt { } dR) 

(35) 

Making use of the classical equation of motion, Eq. (ITrll) , 
Eq. ([33)) becomes 

/ a 2 * /a 2 * a' 2t a _ n a* .„„, 

2m ac 2 2m dd 2 2 2 dt K ' 



2m ac 2 2m ad 2 2 
.dip(c,d,rj) 



1 <9 2 1 3' m 2 2 m 2 2 



C?T7 



ip(c,d,ri) 
(37) 



where 



1 



dt'/ = - / dt' i 



2GM Q 2 
i2 + R 2 



^ (38) 



and 
2 



w (77) 



A 



i 



ml- R s /R + Q 2 /R 2 1 - R s /R + Q 2 /R 2 

(39) 

where w 2 = K/m. Here we have chosen to set ij(t = 0) = 
0. 

To solve Eq. ([3"7]) . we make the ansatz that the wave- 
function is separable, hence we can write the total wave- 
function as 



*(c, d,rj) = ip(c,n)ip(d,ri). 



(40) 



Therefore each of the mode equations satisfies the time- 
dependent harmonic oscillator equation separately. The 
solution for each mode in Eq. ([3"T[) is given by (see 
Ref. [13) 



7Tp 



1/4 



exp 



and 



iJ S , R (d, V ) = e^> — 



1/4 



exp 



vm I p^ 
2 [p 



vm I 
2 \p 



(41) 



7 U 



(42) 

where p v = dp/drj is the derivative of the function p(rj) 
with respect to 77, and p is the real solution to the non- 
linear auxiliary equation 



p vv +^ 2 {v)P = "4 

and where a is the phase, which is given by 

1 f v drf 



a{rj) = -- 



(43) 



(44) 



B. Radiation Only 

In this section we only consider the particles created 
during the collapse of the domain wall. To describe this 
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we use only Eqs. ([27]) and (|28|). In this case the Func- 
tional Schrodinger equation becomes 

~ 2m dc 2 2m dcP + 2 + 2 dt ' [ ' 

Dividing through by / and writing in the standard form 
yields 



1 a 2 * 1 a 2 * m 2t m ,,.,, 2t 



where 



and 



i - + Q 2 /i? 2 ' 



.9* 

(46) 



(47) 



(48) 



where uiq is defined as before. Again separating the wave- 
function as in Eq. (|40p . the solution to Eq. (|4"6")l is the 
same as Eqs. (|4T|) and (|42]) except with the replacement 
r\ now given by 77 (Eq. (1471) ). i.e. 



1/4 



exp 



and 



1/4 



exp 



™ / Pfi_ 
2 U 



2 U 



(49) 



(50) 

Hence each wavefunction evolves with a different time 
parameter. 



C. Taking the Temperature 

For an observer, the complete information about the 
radiation, in the background of the collapsing domain 
wall, is contained in the wavefunction. Consider an ob- 
server with detectors that are designed to register parti- 
cles of different frequencies for the free field $. Such an 
observer will interpret the wave function of a given mode 
at some late time in terms of simple harmonic oscillator 
states. 

For brevity, we consider the following analysis using 
the notation for the entire system. However, the analysis 
for the radiation alone is the same. 

The wave functions in Eqs. (glj) and (gSJ (Eqs. ([4"9]) 
and ([50])) can be decomposed into suitably chosen vac- 
uum basis wave functions, denoted by {</> n } at the final 
frequency u>. The number of quanta in eigenmodes c and 
d can be evaluated by decomposing the wavefunction in 
terms of the states {ip n } and by evaluating the occu- 
pation number of that mode. We decompose the wave- 
function at a time t > tf in terms of a simple harmonic 



oscillator basis at t = 0. The decomposition is 

$(c,d,t) =J2 C n( t ) C m{t)<Pn{c)(p m {d) (51) 



where 



= / dc<fn(c)ip(c,t) and c m = / ddip^ n (d)tl)(d,t) 



(52) 

which arc the overlaps of a Gaussian with the t = 
simple harmonic oscillator basis functions. It was shown 
in Ref. [l6| that the expectation values of the occupation 
number at given eigenfrequency cD (i.e. in the eigenmodes 
c and d) by the time t>tf, are then given by 



N c (t,Q) = ^2n\c n \ 2 and N d (t,Q) = ^ 



m \c m \ 



(53) 

The total occupation number is then the sum of the in- 
dividual occupation numbers, N = N c + Nd- Since the 
initial state is chosen to be the vacuum state, the num- 
ber of quanta in eigenmode measured by the occupation 
number is a consequence of only the collapse, since the 
observer will initially measure zero quanta. 

As stated previously, the observer will interpret the 
wavefunction of a given mode in terms of simple har- 
monic oscillator states. Thus the basis functions (f> n are 
chosen to be simple harmonic oscillator basis states at a 
frequency ui 



(f>n (c) 



\/2 n nl 



(54) 



where H n are the Hcrmitc polynomials. For brevity we 
have only written the simple harmonic oscillator states 
for the c modes, a similar expression holds true for the 
d modes as well. Therefore using Eq. (l52l a nd Eq. ([54]) 
we can write the inner product as in Ref. [ill. [l2l [Tq , for 
n = odd c n = and for n = even 




(wp 2 )V4 



where 



P P 



(55) 



(56) 



Similar expressions hold true for the d modes as well. 

Substituting Eq. (|55|) and Eq. (|56j) and performing the 
sum over n = even in Eq. (|53j) we then have for the c 
modes 



N e (t,w) 



V2 

and for the d modes 

„2 



N d {t,u) 



7f 



i 

top 2 



ujp* 



Pn_ 

Hip 



Pn_ 

ujp 



(57) 



(58) 
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Therefore the total occupation number is given as 



We note that we can rewrite the invariant as 



N = N C + N d = 



2uip 2 



1 - 



En 

Qp 



(59) 



by virtue of Eqs. (|57|) and ([58]) . By fitting the occupation 
number to that of the Planck distribution, 



N P 



1 



; /3w _ i 



(60) 



the occupations at eigenmodes c and d can then be used 
to find the temperature of the radiation. From Eqs. (|4"Tj) 
and (gSJ) and Eqs. (g5J) and (J5DJ) we see that the occu- 
pation number for each of the systems will be different, 
since both rj and fj are different. 



IV. ENTROPY FUNCTION 

Here we develop the entropy function for studying the 
time evolution of the collapsing domain wall. Again for 
brevity we use the notation for the entire system, where 
again the analysis is the same for the radiation only sys- 
tem. 



Partition Function 



To study the entropy of the system, we will first de- 
velop the partition function for the system. Following the 
procedure used in Ref. for the Liouville-von Neumann 
approach, we can write the partition function as 



Z = Tr [e-? 1 ] 



(61) 



where / is any operator which satisfies the Liouville-von 
Neumann equation Eq. (J2|) and [3 is a free parameter. 
From Ref. fl7| . we can write the invariant operator I 

as 



1 



+ (n c p - mp„c) 2 +\ - + {TTdP - mp„d) 2 

p V p 



We can therefore write the partition function as 
,1 



(62) 



Z = Tr exp 



ft' 



c 1 { \ 2 
+ {TT c p - mp„c) 

P 



- + {TT d p - mp v d) 
P 



(63) 



1= ^ 



1 

V2 



1/4 



- i (ir c p - mp v c) 



1/4 



+ i (n c p - mpr,c) 



1/4 



1/4 



i (ir d p - mp v d) 



where 



and 



+ i (ir d p - mp v d) 
=n(t) + m(t) + 1 

n(t) = a t (t)a(t), and m(t) = a\t)d{t) 



(64) 



(65) 



-(*) - T2 



1/4 



+ i (n c p - mpnc) 



ait) 



V2 



1/4 



+ i {-KdP - mpr,d) 



(66) 



Here n(t) is the time-dependent number of states for the 
c modes, while m(t) is the time-dependent number of 
states for the d modes. At a particular time t, one has in 
Fock space 

n(t) In, m, t) — n\n, m, t), and m(t)\n,m,t) = m\n,m,t). 

(67) 

Thus, in this space we can then write the partition func- 
tion as 



Z 



Tr exp [— /3u;o (n + m + 1)] 
1 



4 sinh' 



(68) 



Here we note that Eq. (|68|) takes the form of a two- 
dimensional harmonic oscillator. 

In Ref. [HI, [12, EBl one can define the occupation num- 
ber for a frequency w, Eq. (|59|) . Then by fitting the num- 
ber of particles created as the usual Planck distribution 
Eq. (|60[) . one can then in principle fit the temperature of 
the radiation. Here, we then choose to define j3 as 



P = 



dhx{\ + 1/N) 
dui 



(69) 



This implies that all of the time-dependence of the system 
is encoded into the temperature of the system. 

We can see that Eq. (|68|) is just the standard entropy 
for a time-independent harmonic oscillator; however, the 
temperature here is time-dependent. Therefore the en- 
tropy is also time-dependent. 



7 



B. Entropy Analysis 



Using Eq. (|68[) and the thermodynamic definition of 
entropy, we can then write the entropy of the system as 



S = -ln((l-e-^ ) 2 ) +2(3- 



(70) 



Therefore, this is again just the entropy of the usual two- 
dimensional time-independent harmonic oscillator. From 
Eq. ([5]) it follows that the temperature is time-dependent. 

First we consider the entropy of the entire system. In 
Figure [T] we plot the entropy of the entire system as a 
function of dimcnsionless time tf-/R+. Figure Q] shows 
that the system starts with an initial entropy of zero. 
This is expected since initially there is only one degree 
of freedom, meaning that S = ln(l) = 0. Here we have 
normalized the initial entropy of the shell to be zero. To 
justify this normalization, consider a solar mass black 
hole. Under the usual Bckcnstcin-Hawking entropy, the 
order of magnitude estimate of the entropy of a solar 
mass black hole is Sbh ~ 10 75 . Now consider that the 
shell is actually made up of protons. The initial entropy 
of the shell then is approximately Ss,o ~ 10 57 . Compar- 
ing the entropy of the final black hole versus the initial 
entropy of the shell, the entropy of the final black hole 
is much much greater than that of the initial entropy of 
the shell, thus the initial entropy of the shell only con- 
tributes a negligible amount of entropy to the entropy of 
the final black hole. Thus our normalization of the ini- 
tial entropy of the shell to zero is justified. As tf^/R + 
increases, initially the entropy increases rapidly, then set- 
tles down to increase approximately linearly. Due to the 
linear increase, we see that as tf-/R+ goes to infinity, the 
entropy will then diverge. This is again expected since 
as the asymptotic time goes to infinity, the number of 
particles that arc produced diverges (see Ref. [l6| ) . This 
is a consequence of the fact that we keep the background 
fixed (i.e. i?+ is a constant). In reality, R+ should de- 
crease over time since the radiation is taking away mass 
and energy from the system. Therefore as t/_/i? + goes 
to infinity, the entropy of the entire system as measured 
by the asymptotic observer diverges as R — > i? + . 

Now we consider just the radiation which is induced 
during the collapse. In Figure [5] we then plot the en- 
tropy as a function of t-/R + . Figure [5] shows initially 
the entropy of the system is zero. Again, this is expected 
since initially the domain wall is in vacuum, meaning that 
initially there is no radiation being induced. Therefore 
the only degree of freedom is that of the domain wall, 
which then gives that the initial entropy of the radia- 
tion must be zero. As increases, initially there 
is a rapid increase in the entropy, but again, the entropy 
then increases linearly as tf-/R+ increases. As in the 
case of the entire system, as tf-/R + goes to infinity, the 
entropy of the radiation also diverges. This is expected 
since the induced radiation diverges as R — > i? + ; hence as 
the domain wall approaches the horizon, the occupation 
number for the induced radiation diverges. 




t/_/R+ 



FIG. 1: We plot the entropy of the entire system as a func- 
tion of tf-/R+. Except for the initial increase, the entropy 
increases approximately linearly as tf-/R+ increases. 



In Figure [3] we plot the entropy as a function of asymp- 
totic observer time t of both the entire system and the 
induced radiation during the time of collapse. Fig- 
ure [3] shows that except the initial increase in the en- 
tropy, for later tf-/R+ (approximately tf-/R+ = 7.5), 
the slopes of the entropy versus time are approximately 
equal. Therefore, one can expect that the entropy of the 
domain wall is approximately constant for late times. 

What is of interest is the entropy of the collapsing 
charged domain wall, since this will form a black hole. 
To find the entropy of the charged domain wall, we can 
take the entropy of the entire system and subtract off the 
entropy of the induced radiation (since these are the only 
relevant objects which contribute to the entropy). The 
result is then given in Figure |4j We see some interesting 
features of the time-dependent entropy in Figure |4] Here 
we discuss these features; first we will make some general 
comments on the overall behavior of the time-dependent 
entropy, followed by discussions of late-time (values of 
tf-/R+ > 8) and the early-time (values of t-/R+ < 8) 



S(t) 
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FIG. 2: We plot the entropy of the induced radiation created 
during the collapse as a function of tf-/R+. Except for the 
initial increase, the entropy increases approximately linearly 
as tf-/R + increases. 
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FIG. 3: We plot the entropy as a function of tf-/R+ for both 
the entire system and the induced radiation during the time 
of collapse. At the time tf-/R+ ~ 7.5, the slope of the two 
lines is approximately equal. 

behaviors of the time-dependent entropy, respectively. 

Figure 2] shows that initially the entropy of the charged 
domain wall is zero. As stated above, this is expected 
since initially there is only one degree of freedom. As 
tf-/R+ increases, the entropy of the domain wall rapidly 
increases. However, for late times (tf-/R+ > 8), the 
entropy of the charged domain wall goes to a constant. 
From the discussion above, this is expected since the late- 
time entropies for the entire system and for the radiation 
are approximately parallel. This feature is somewhat 
artificial though: the entropy is constant since we are 
assuming that the mass is approximately the Hamilto- 
nian of the system, which is a constant of motion (see 
Ref. [H[). This means that since we are holding the 
mass of the charged domain wall constant, we need to 
keep adding energy to the system to counteract the loss 
of mass (energy) from the radiation. Therefore one can 
expect that the entropy of the charged domain wall must 
be a constant for late times. 

In reality, radiation takes mass away from the system, 
so the entropy of the charged domain wall will go to zero 
as R + goes to zero. Hence if the charged domain wall 
doesn't start off with enough mass, it could evaporate 
before ever creating the black hole. This means that after 
the charged domain wall disappears, all the entropy will 
go into the entropy of the radiation, since all the mass 
has dissipated. 

For large values of tf-/R+ (tf_/R + > 8) the entropy 
of the charged domain wall is constant. From Figure 0] 
we see that our numerical value for the entropy of the 
charged domain wall for this region is 

S«2(i?+//_) 2 . 

Comparing with Eq. |T]), we can view this discrepancy 
as a shift in the horizon radius In order to get the 
theoretical value for the entropy, Eq. ([j}, we see that 
we would require i?+ — > 1.25i? s . This is an understand- 
able numerical error, which implies that our numerical 
solution is of the same order as the Hawking-Bekenstein 



Entropy of Charged Shell 

S(t) 




FIG. 4: We plot the entropy of the charged domain wall as a 
function of tf-/R+. Here we see for tf-/R+ > 8 the entropy 
is constant in time with a value of S ~ 2(i? + /_/L) 2 . 

entropy This is in better agreement than the result for 
the Schwarzschild case, see Ref. 

For small values of tf-/R + (tf-/R+ < 8), the en- 
tropy of the charged domain wall seems to stop increas- 
ing around ps 7.5. This means that the entire 
change in entropy of the domain wall occurs in the region 
< tf-/R+ ps 7.5. At first this value for tf-/R+ seems 
arbitrary, however there is some physical significance to 
this time. 

First, we notice that the volume of the charged domain 
wall becomes approximately constant at tf-/R+ ss 7.5. 
To illustrate this we will first require Rq = nR + , where 
n g S (S being the real numbers). From Eq. (f2"2"j) we can 
then write the position of the domain wall as 

R{t) =R+(l + (n - l)e- tf -/ R +^j . 

For illustration purposes we choose the value n = 10, 
which at time = 7.5 corresponds to a distance 

of R ps 1.005i?+ with velocity R = -f ps -0.005 (the 
numerical values for the position and velocity will be 
even less for smaller values of n). Hence the value 
of tf-/R+ ps 7.5 corresponds to the region where the 
charged domain wall is very close to the horizon and the 
velocity of the charged domain wall, with respect to the 
asymptotic, goes approximately to zero. Therefore, as far 
as the asymptotic observer is concerned, the dynamical 
process of the charged domain wall collapsing is over and 
the domain wall has come to rest, meaning that there is 
now approximately a constant volume. From this point 
on, classically it takes an infinite amount of time for the 
domain wall to reach the horizon radius, hence travel the 
distance R ps 0.005i?+, see Section ITT1 

Second, we can show that the entire system and the 
induced radiation come into thermal equilibrium at this 
time. In Figure [5] we plot /3 versus tf-/R + for the en- 
tire system (continuous curve) and the induced radia- 
tion (dashed curve) . Figure [5] shows that for the time 
tf-/R+ ss 7.5 the values of the two /?'s become approxi- 
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almost linearly, except for the bump or oscillation occur- 
ring during the time span 5 < tf_/R + < 10. This oscilla- 
tion may be attributed to several different circumstances. 
First, the oscillation may be caused by the non-thermal 
property of the radiation (see Ref. [16[ ) ■ Secondly, this 
oscillation may be a manifestation of the error associated 
with the numerical calculations. Lastly, the oscillation 
may be an artifact of expanding the calculations beyond 
the region of validity, since we are using the near horizon 
approximation. Hence for values large compared to R + , 
we cannot completely trust our result. 



FIG. 5: We plot ft of the entire system and the induced radi- 
ation as a function of tf-/R+. Here we see for tf-/R+ ~ 7.5 
the two values of /3 become approximately equal, meaning the 
two systems are in thermal equilibrium. 

mately equal, meaning that the entire system and the in- 
duced radiation are now at the same temperature. There- 
fore the system is now in thermal equilibrium, meaning 
that there is no more change in entropy of the charged 
domain wall as tf-/R + increases. Further more, the fluc- 
tuations (departure from thermalit y) i n /3 become smaller 
at this time, as discussed in Refs. jl6j . 

Finally we can evaluate the the chemical potential for 
both the entire system and for the induced radiation. 
From definition we can write the chemical potential as 



In Figure |5] we plot the chemical potential for both the 
entire system and for the induced radiation. We can see 
that as tf-/R+ increases the chemical potential of the 
entire system and the induced radiation goes to zero. 
This means that the dispersion of particles goes to zero 
and the system goes into equilibrium. 




FIG. 6: We plot /i versus tf-/R+. The solid line corresponds 
to the entire system while the dashed line corresponds to the 
induced radiation only. Here we see that as tf-/R+ increases, 
the chemical potential for each goes to zero. 

During the dynamical process, the entropy increases 



V. CONCLUSION 

We studied the time evolution of the entropy of a col- 
lapsing spherically symmetric charged domain wall (rep- 
resenting a shell of charged matter) using the Func- 
tional Schrodinger and Louisville-von Neumann for- 
malisms from the point of view of an asymptotic observer. 
This was done by investigating the change in entropy for 
the entire system (the charged domain wall and the in- 
duced radiation) and that of only the radiation induced 
during the collapse semi-classically, i.e. the background 
of the collapsing charged domain wall was treated classi- 
cally while the radiation was treated quantum mechani- 
cally. As stated in Section HH this was done for the case 
(GM) 2 > Q 2 only. The reason for this is that in the 
case Q 2 > (GM) 2 Eq. (|22[) contains an imaginary term 
in the exponential term. Therefore the solution contains 
sinusoidal terms, where the charged domain wall will os- 
cillate at some radius R outside the horizon radius. Thus, 
the horizon is never formed and a black hole never gets 
created. 

To summarize, we have arrived at the following re- 
sults. From Figure[T]we see that dS/dt > for the entire 
system, which is in agreement with the second law of 
Thermodynamics. Second, from Figure [5J we see that 
initially the entropy of the radiation increases at a lower 
rate than that of the entropy of the entire system; how- 
ever, for larger values of tf-/R+ the increases of each 
are approximately the same (as seen in Figure [3]). The 
difference in the initial rates is due to the fact that each 
system evolves with a different time parameter (77 and fj). 
Finally, from Figure 2] wc demonstrated that the entropy 
of the charged domain wall is constant for late times. 
As discussed in Section IIVB1 this is expected since we 
are holding the mass of the charged domain wall con- 
stant, thus we must keep adding energy to the system to 
maintain the mass of the charged domain wall. In a com- 
pletely realistic model, this would not be the case. The 
entropy would be expected to decrease as R+ decrease, 
since R+ depends on the mass of the charged domain 
wall. So as the mass decreases the horizon radius would 
also decrease, meaning that the entropy of the charged 
domain wall would go to zero as R+ goes to zero. Here 
we demonstrated that the late-time entropy (large values 
of tf-/R+) of the charged domain wall does in fact go 
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to a constant as predicted by Bekcnstein, see Refs. [3|, |4|. 
We found that the late time entropy of the charged do- 
main wall is S « 2(i? + //_) 2 . By comparing Eq. (fTJ) with 
our result, we see that our result is consistent with the 
accepted value for the entropy. The entropy for a col- 
lapsing massive shell was also derived in Ref . [l9| , where 
the authors used a mini-superspace effective action for 
the area and mass aspect of the shell. The authors show 
analytically that the entropy of the shell is on the order 
of the Bekcnstcin-Hawking entropy. 

Figure |4] displays some other interesting features. First 
we see that all the change in entropy occurs during small 
values of tf- / R+ , with the increase in entropy ending at 
a time tf-/R + rj 7.5. As discussed in Section HVB| this 
time is not an arbitrary time. To see this, first, for differ- 
ent values i?o, this time corresponds to when the charged 
domain wall is very close to the horizon and the corre- 
sponding velocity is very close to zero. From this point 
on the charged domain wall takes an infinite amount of 
time to traverse the remaining distance to the horizon. 
Thus, as far as the asymptotic observer is concerned, 
the dynamics of the charged domain wall has finished by 
the time tf-/R + « 7.5, hence the observer won't sec 



any further change in the entropy since the volume of 
the charged shell is approximately constant. Second, the 
"inverse temperature" (/3) for both the entire system and 
the induced radiation become equal at this time. This 
means that the two systems come into thermal equilib- 
rium, as seen in Figure [5] Furthermore, the fluctuations 
in the temperature become smaller by this time, as dis- 
cussed in Rcf. [l6| . However, these oscillations may still 
be substantial and be the cause of the oscillation present 
from 5 < tf-/R+ < 10 (see Figured]). Finally from 
Figure [6] we see that the chemical potential for both the 
entire system and the induced radiation go to zero near 
this time. Again, this signifies that the dispersion of the 
induced particles goes to zero and the system arrives at 
equilibrium. 
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